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Abstract. Surprisingly, differentiable functions are able to oscillate arbitrarily faster than their 
highest Fourier component would suggest. The phenomenon is called superoscillation. Recently, 
a practical method for calculating superoscillatory functions was presented and it was shown that 
superoscillatory quantum mechanical wave functions should exhibit a number of counter-intuitive 
physical effects. Following up on this work, we here present more general methods which allow the 
calculation of superoscillatory wave functions with custom-designed physical properties. We give 
concrete examples and we prove results about the limits to superoscillatory behavior. We also give 
a simple and intuitive new explanation for the exponential computational cost of superoscillations. 
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1. Introduction 

It used to be believed that a function could not oscillate much faster than its highest 
Fourier component. Aharonov, Berry and others showed that this is not the case by giving 
explicit counter-examples which they named superoscillatory functions, see, e.g., [I]- 111- In 
fact, there are functions which on arbitrarily long stretches oscillate arbitrarily faster than 
their highest frequency Fourier component, see [Sj- In other words, the presence of localized 
fast oscillations in a continuous function need not be visible at all in the function's global 
Fourier transform. In a function's global Fourier transform, contributions from regions of 
fast oscillations can be cancelled perfectly by contributions from regions where the wave 
function is oscillating slowly. 

In the context of quantum theory, wave functions that superoscillate are able to cause 
a number of counter-intuitive effects. Some of these may be of conceptual significance in 
quantum gravity, see H] • But effects of superoscillations also enter in the low energy realm 
of nonrelativistic quantum mechanics. Among such potentially observable low-energy effects 
is the counter-intuitive phenomenon that particles with superoscillatory wave functions can 
be made to accelerate when passing through a neutral slit: 

Consider a particle which possesses a bounded momentum range, i.e., its momentum wave 
function vanishes for momenta that are larger than some p max . As will be explained below, 
we can arrange that in a certain region in space the particle's wave function superoscillates, 
i.e. that it oscillates with a much shorter wavelength than h/p max . Now let the wave function 
be incident onto a screen with a single slit in such a way that it is the superoscillatory part of 
the wave function which passes through the slit. Upon emerging from the slit the particle's 
wave function will then oscillate rapidly where the slit is and will be zero elsewhere. The 
very short wavelengths of the emerging wave function will be visible in its global Fourier 
transform. This is because the contributions to the global Fourier transform which come 
from the fast oscillations in the slit interval are no longer cancelled by contributions from 
outside the slit interval. Therefore, the particle will have gained momentum merely by 
passing the slit. The momentum gain is determined by the shortness of the wave-length of 
the superoscillations and, as explained below, there is no limit, in principle, to how short 
that wavelength can be made. 

In order to facilitate the design of experiments that can realize the effects of superoscilla- 
tory wave functions it is desirable to possess methods for explicitly calculating superoscilla- 
tory wave functions with predetermined properties. In particular, one may wish to calculate 
those low-momentum but superoscillatory wave functions which after passing through the 
slit yield wave functions with a predetermined arbitrarily large momentum and a momentum 
uncertainty that is as small as is allowed by the uncertainty relation. Our aim here is to 
develop methods that allow us to solve this and other problems. 

Our starting point will be the method for calculating superoscillatory wave functions which 
was developed in using results of and 0. This method allows the construction of wave 
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functions of arbitrarily low fixed frequency content that pass through an arbitrary finite 
number of pre-specified points. Figures ^ and 121 show an example. 




Figure 1 . Example of a superoscillation created by requiring the wave func- 
tion to pass through certain points. A cosine function of the minimum wave- 
length/maximum frequency is shown for comparison. 




Figure 2. A zoomed-out figure^ Notice that, as is typical, the amplitudes in 
the superoscillating region are far smaller (here even unnoticable) than those 
on either side. 

Our aim is to develop more general methods for designing superoscillatory functions with 
generic pre-specified properties. We will also ask what the in-principle limits are for the 
construction of superoscillatory wave functions. 

2. Self-acceleration through single slit 



In order to motivate and formalize the mathematical problem that we will address, let us 
consider the illustrative example of particles that self-accelerate when passing through a slit. 



4 



2.1. Notation. We will denote a wave function tp and its Fourier transform by the same 
letter, since they correspond to the same vector in the Hilbert space of states. 

1 f°° ipx 

ip(p) = —:= / ip(x) e~~ dx (1) 



In the text, whenever necessary, we will write ip(x) or ip{p) to indicate if the position or the 
momentum wave function is meant. We will often consider particles whose momentum is 
bounded by a finite value p max : 



Pmax 



if)(x) = / i>(p) e k dx (2) 

V2irh i- Pmax 

Borrowing terminology from communication engineering and sampling theory, see e.g. |lUj . 
we will speak of such a function ip(x) as having bandwidth p max , as being band-limited, or 
in this case as being momentum-limited. It will be convenient to define the sine function as: 

sine (x) := < 3 
\l, if x = 

Notice that definitions of the sine function elsewhere may include a factor of tt. 

2.2. Gedanken experiment. Let us consider a particle in two dimensions which travels 
along the Xi-direction towards a screen which is parallel to the ^-direction. Assume the 
particle passes through a slit with a^-coordinate interval [—§,§] in the screen. Henceforth, 
we will assume that the incident particle's momentum parallel to the screen, p 2 , has a finite 
bound p 2max : 

^(pi,P2,t) =0 if p 2 ^ [-P2 max ,P 2max ] (4) 
Our aim is to compare the particle's momentum parallel to the screen before and after the 
particle passes the slit. For simplicity, we will suppress the variables x\,p\ and t. From now 
on, X2 is renamed x and p<± is renamed p. We denote the incident wave function just before 
passing through the slit by if)(x) and we denote the wave function which emerges from the 
slit by *f?(x). The state \^f) is of course given by projection and renormalization 

II p s W) II 

where P s projects onto the slit: 

P s ^{x) := < . (6) 

|^U, otherwise 

Similarly, we define P& as the projector onto a finite momentum interval: 

D U \ JV(p), if \P\ <Pmax 

PbnP) ■= < n ^ , (7) 
[(J, otherwise 

While the incident wave function is momentum limited, Pb\ip) — the emerging wave 
function is position limited, obeying P s \^) = \^)- 
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As was shown in E] , it is always possible to find incident wave functions ip(x) which 
obey the momentum bound p max and which at any finite number of points in the slit interval 
take arbitrarily prescribed amplitudes (we will reproduce this result as a special case below). 
We will be able to arrange, therefore, that the wave function ip(x) takes for example the 
alternating values (— l) n at an arbitrarily large number of points in the slit interval - which 
enforces superoscillations. These ijj(x) will be different iable and square integrable. Then, if 
the particle passes the slit, only the superoscillating stretch of the wave function emerges 
from the slit. Renormalized, we denote it *f>(x). The Fourier transform of *f?(x) will 
show the presence of small wavelengths, implying that the particle emerges from the slit 
accelerated to a momentum beyond p max - 

2.3. Template Functions. As already mentioned, the results of [SlIE] showed that functions 
of fixed bandwidth can always be found which at arbitrarily but finitely many points possess 
predetermined amplitudes. Therefore, the width or narrowness of the slit does not limit 
how short the wavelength of the superoscillations can be. As a consequence, there is no 
slit-dependent limit to the amount of self-acceleration that can be achieved in this way. 

This leads us to ask more generally whether the process of self-acceleration can be designed 
virtually at will: is it always possible to construct incident wave functions ip(x) of fixed mo- 
mentum limit p m ax which on the slit interval [ — -§] match any arbitrarily-chosen template 
function, say $(x)? This is of interest because, if true, we can optimize the predictability of 
the self-acceleration. To this end, we would choose the template function to be a wave 
function with a fixed arbitrarily large momentum expectation p whose momentum uncer- 
tainty Ap is as small as allowed by the uncertainty relation. If the incident superoscillatory 
wave function matched this template function in the slit interval (up to normalization), then 
the wave function would merge from the slit with the chosen momentum expectation p and 
lowest possible momentum uncertainty, Ap, for the given width of the slit. For later use, let 
us calculate these ideal template functions <&(x). 

2.3.1. Ideal Template Functions. Our aim is to find 'ideal' template functions $(x) defined 
on the slit interval [ — §, §] which minimize the momentum uncertainty Ap, possess a prede- 
termined momentum expectation ($|p|$) = p and are normalized ($|$) = 1. To this end, 
we need to solve the constrained variational problem with the functional 



where /ii and [1% are Lagrange multipliers. Note that minimizing Ap is equivalent to mini- 
mizing ($|p 2 |$) because p is fixed. Hence the Euler-Lagrange equation in position space is: 



Since any wave function that emerges from the slit vanishes at the slit boundaries, we require 
$(±L/2) = 0. The solution is unique up to a phase: 






(9) 




(10) 
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Its uncertainties are Ax = L((tt 2 — 6) /12-k 2 ) 1 / 2 m 0.18L and Ap = We have AxAp « 
0.57ft which is a little larger than what the uncertainty relation allows because our problem 
requires template functions <f>(x) to be zero outside the slit interval. 

2.3.2. Superos dilatory Wave Functions Cannot Match Arbitrary Templates. Let us now 
come back to the question whether it is generally possible to find an incident wavefunction 
ip(x) which obeys a momentum bound p max while in the interval [— -§] agreeing completely 
with an arbitrarily chosen template function <&(x), such as the function just calculated 
in Sec J2.3.ll Strictly speaking, the answer is no: 

As is easily verified, all bandlimited functions are entire functions. In particular, any 
momentum-limited incident wave function ip(x) is entire and it is, therefore, everywhere 
differentiable. Now choose, for example, a template function which is not differentiable 
at some point in the interval [— -|, |f] . Thus, there cannot exist a momentum- limited incident 
wave function which obeys ip(x) = for all x G [— |f , ^] . Nevertheless, a slightly weaker 
proposition does hold. 

2.3.3. Convergence Towards Arbitrary Template Functions. Let <&(x) be a continuous and 
square integrable template function. Let us ask whether one can always find a sequence of 
wave functions iPn{x) of fixed momentum bound p max which behave with more and more 
precision like over the region of the slit. To be precise, is it possible to find a se- 
quence of momentum-limited incident waves \iPn) whose emerging wave functions \^n) have 
asymptotically vanishing L 2 -distance HI^jv) — 1^)11 to an arbitrary template state |$)? This 
is indeed the case. 

To see this, consider in the quantum mechanical Hilbert space of states TC with scalar 
product 



and the following three subspaces: 





/■oo 

/ « 

' — oo 


'(x)^2{x) dx 


n s 


:=P S 


n 


TCb 


■=p b 


n 


7~Lsb 


■=Ps 





(11) 

(12) 
(13) 
(14) 

That is, H s is the subspace of states with position-limitation to the slit, Tib is the subspace 
of states with fixed momentum-limitation p max , and 7i s b is the subspace of states obtained 
after passing the momentum-limited wave functions through the slit. 

Proposition 1. H s b is dense inH s , i.e.: 

V |$) G H s ,e > 3 |tf> G H sb : |||tf> - \$)\\ < e. (15) 

Proof. If |$) = 0, then take \^) = 0. For |$) ^ we must show that: 

$ |$) G Hs\{0} : (*|$) = V |*) G H sb (16) 



7 



Since |$) is position-limited, this is equivalent to showing that: 

$ |$> e H s \{0} ■■ (V#> = o V |^> e n b (17) 

Assume, for a contradiction, that: 

3 |$) G H s \{0} : (V>|$) = V |V) e % (18) 

This implies that |$) _L Thus, $(p) = on [— J9 max ,p max ]. But, since $(p) is entire and 
zero over a finite interval, = everywhere on M, i.e., |$) = 0. This is a contradiction. 
Therefore, is dense in 7i s . □ 

While this result proves the existence of bandlimited functions that are arbitrarily close 
in the L 2 topology to any template function within the window of the slit, the result does 
not provide explicit methods for constructing such bandlimited functions. 

3. Constructive Method for General Linear Constraints 

We now focus on practical methods for calculating superoscillatory wave functions that 
approximate template functions in the slit interval. We begin with the method for con- 
structing superoscillatory functions presented in [Hj. This method allows one to specify that 
the to-be-found superoscillatory function takes arbitrarily chosen amplitudes a k at any finite 
number N of arbitrarily chosen points x k : 

^ u \x k )=a k for k = 1,..., N (19) 

The superscript ^ is to indicate that the function will generally be unnormalized. EqJT9l 
specifies a function which possesses a superoscillating stretch. For example, we may choose 
the x k spaced closer than h/p max and the amplitudes alternating, e.g. a k = (— l) k . The 
normalized wave function ip(x) = ^ u \x)/\\ip^\\ then possesses superoscillations that are 
as rapid as those of ip 1 ^ but with a renormalized amplitude. Thus, in order to obtain the 
ip(x) with the most pronounced superoscillations, i.e. the superoscillations of largest possible 
amplitude, one needs to find that function ip^(x) whose L 2 norm ll"?/^!! is minimal. The 
method of jHj solves this optimization problem. 

We now generalize the method of To this end, we begin by rewriting the requirement 
that ip( u \x) be bandlimited by p max and pass through the points {(x k , Qfc)}fcLi> namely (|TU)l. 
in momentum space as 

1 /-Pmax 

-7= / e i **i/>M(p)dp = a k . (20) 

Our aim is to obtain a method for constructing superoscillatory wave functions which not only 
pass through predetermined points but which obey also more generic types of constraints. 
To this end, let us allow constraints on the function vjj^ which are of the general linear form: 



V2nE 



Pmax 



Xl{p)^ u \p)dp \/ke{l,...,N} (21) 



Pmax 



Here, the \k are arbitrary linearly independent differentiable functions. By choosing these, 
we will be able to prescribe for the superoscillatory wave function not only amplitudes but 



8 



also arbitrary derivatives, integrals and any other linear constraint. In order to obtain the 
most pronounced superoscillations in the normalized function ip we minimize the norm of 
ip( u \ subject to the constraints in Eql21l The to-be-optimized functional with Lagrange 
multipliers A& reads 



'max 



C = / ^\p)^ u \p) dp-J2 -7= / Xl(p)^ iu) (p) dp + c.c, (22) 
leading to the Euler-Lagrange equation: 

1 N 

V> ( %) = ^=£A fe x fe (p) ( 23 ) 



Recall that ip^ is zero outside the interval [— p ma , x ,p m & x } by assumption. Thus, using (f2~3J) 
in (j2U), 



N 

a fc = ^T fcr A r , (24) 



r=l 



where the Hermit ian matrix T is defined by: 



T kr := — — I xl(p) Xr(p) dp (25) 
27m /__ 

" /'max 

As we will show below, T is invertible. Thus, A = T~ l a, i.e.: 

N 



r=l 

Thus, using the Fourier transform of the constraint function 

Xfc(z) : = ^ / X*(P) e ^ rf P ( 27 ) 

we obtain from (J23J) that the desired superoscillatory (still unnormalized) incident wave 
function in position space is given by: 



1 N 

^ u) (x) = ^=J2^Xk(x) (28) 



k=l 



3.1. Existence of the Solution. It remains to be shown that T is indeed invertible. To 
see this, let u be an arbitrary vector. Then: 

N 



Since the Xk are linearly independent the integrand is positive. Therefore, T is positive 
definite and hence invertible. 



As was shown in [6 , one cost of superoscillations is that requiring more or faster super- 
oscillations makes the matrix T increasingly difficult to invert numerically, as its smallest 
and largest eigenvalues differ by growing orders of magnitude. The condition number was 
found to increase exponentially with the number of superoscillations. 

We here only remark that, in the sense of computational complexity, this makes it compu- 
tationally hard to calculate superoscillations. Interestingly, this also means that any quantum 
effect that naturally produces functions with arbitrarily large superoscillatory stretches con- 
stitutes an example of an exponential speed-up in the sense of quantum computing. Physical 
occurrences of superoscillations, e.g. in the context of evanescent waves, have been discussed 
e.g. in Also, for example, (rather speculatively) the possibility of an unbounded pro- 

duction of superoscillations has been discussed in the context of the transplanckian problem 
of black holes in [HI EES- 

Here, we will focus on a more immediate cost of superoscillations, namely the need for 
an increasingly large dynamical range: a function's superoscillations are generally of low 
amplitude when compared to the function's amplitudes to the left and right of its super- 
oscillatory stretch. To be precise, it was shown in jB] that the L 2 norm of the function 
increases polynomially with the frequency of the prescribed superoscillations, for fixed pre- 
scribed superoscillating amplitudes. In particular, it was also shown that the norm increases 
exponentially with the number of imposed superoscillations. Correspondingly, in normalized 
wave functions the amplitudes of superoscillations decrease exponentially with the number 
of superoscillations. (Of course, if the superoscillating stretch of the particle's wave function 
happens to pass through the slit then its wave function, however small, is re-normalized 
whereby the superoscillating amplitudes will be restored to the amplitudes of the template 
function.) 

By making use of the special properties of prolate functions these scaling results were 
derived for the type of superoscillations produced with method of jH]. In the following two 




fc,r=l 





4. The Cost of Superoscillations 
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subsections we will show more directly the underlying reason for this exponential behavior 
of the norm of superoscillatory functions. Our argument will apply more generally to all 
superoscillatory functions that arise from linear constraints. 

4.1. Derivatives and Norms. If a function ^ u '(x) is bandlimited one would expect that 
there is a bound on its derivatives. Applying the Cauchy-Schwarz inequality, consider 



dx n 



(«), 



x) 




1 dp 



Thus, 



d" 



dx" 1 



(«). 



x. 



< 



h 



Pn 



ixh 



Ml 



(29) 



Thus arbitrarily large derivatives, as they can be produced with superoscillations, are con- 
sistent with a finite fixed bandwidth but we see that the cost must be an increase in the 
norm of the function 

4.2. The Norm of Superoscillating Functions. A precise expression for the norm ||?/>^ || 
of the superoscillatory functions obtained by our method can be derived: 

2 



Ml 



Hence: 



2nh 
l 

2vrfr 

A f rX 

(«)|| 2 - 



pmax 

"pmax 
Pmax 



N 



Y^kXk(p) 



k=l 
N 



dp 



Yl X *kXl(p)Xr(p)K dp 

Pmax f> = \ 



(30) 



(31) 



Note that T _1 is a positive self-adjoint matrix. We now see that for given constraint func- 
tions, Xk, the most norm-expensive superoscillatory functions are obtained if we choose the 
constraint parameters a k such that a is the eigenvector of T -1 with largest eigenvalue. We 
will arrive at those extreme superoscillations also from independent momentum space con- 
siderations in SecEl 

4.3. Adding Successive Constraints. Consider a set of constraints, described by a set 
of functions {xk}k=i an d parameters {a k }^ =1 and suppose that, using our method, the 
momentum-limited wave function which obeys all those constraints and is of minimum norm 
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has been calculated. Let us ask how the norm of the solution to this problem changes if we 
add one additional constraint 

aN+1 = ~7tt^ / XN+i(p)ip(p) dp (32) 

where Xn+i and a^v+i are chosen arbitrarily. Let us denote the solution to the initial problem 
of N constraints by ip^ and let us define: 

Pmax 

Xn+i(p)^Pn(p) dp (33) 

Pmax 



V2irhJ- 

Clearly, if we choose the (N + l)st constraint with cin+i := c then ip^ is also the function of 
minimum norm obeying the N + 1 constraints, i.e. ipN+i = ipN, just as if we had not added 
a new constraint, or as if we had set the (N + l)st Lagrange multiplier to zero: Aat + i = 0. 

Now, let us allow the constraint parameter cin+i to vary away from c. Correspondingly, our 
method will yield a family of functions, ipN+i ?/V), parametrized by ajv+i- We observe 
from ()31|) (letting the sum run to N + 1) that the norm squared of these functions is a 
quadratic (and of course positive) polynomial in a m . Note that its minimum occurs if we 
choose the aAr +1 -value 

C = a N+l = - x Yl T (N+l),r a r ( 34 ) 

I (N+l),(N+l) r ^(N+l) 

because then X^(at+i) ^(tv+i) s° s = ^- Using (|2^j) we see that this choice of aAr+i leads to 
the vanishing of the Lagrange multiplier A^+i = 0, which is what we expected for if we add 
a new constraint that is already satisfied, if) will not change. 

Crucially, we now see that as we tune a^+i away from c, say in order to enforce an 
additional superoscillation twist, the squared norm of the solution increases quadratically. 
Therefore, if we keep adding new generic constraints, say in order to implement more and 
more superoscillations, this will generally increase the norm of the solution by a factor in 
each step. Thus, the norm of the solution will generically scale exponentially with an increase 
in the number of constraints N. 

This finding widely generalizes the result of jH] which applied only to constraints of the 
special form (fTU|) and among them only to those with equidistant spacings of the x^. 



5. Applications to an 'Ideal' Template Function 

In Sec 12.3.11 we asked how the wave function \l/(x) that emerges from the slit would have 
to look in order to describe a particle with an arbitrarily high predetermined momentum 
expectation value p and a momentum uncertainty Ap which is as small as is allowed by the 
uncertainty relation. This 'ideal' template function was given in (jlOj). 

Let us consider the concrete example, h — 1, L — 2ir, p max = 1, and p = 2. If the emerging 
wave function \1/ can be arranged to be equal or close to this template function $, this clearly 
exhibits the phenomenon of self-acceleration because the emerging momentum wave function 
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would be peaked at p ~ p = 2, i.e. well outside the original bandwidth of p 
FigH 



max 1) see 




Figure 3. The 'ideal' template's momentum wave function . Notice 

that it is centred well outside the original bandwidth p max = 1- 



We had shown that exact matching, = $(x), is generally not possible, but we also saw 
that there always exists a sequence of incident waves ipN so that for the emerging waves 
we have ^/jv — > $ in the L 2 topology, which is here the only physically relevant topology. 
Thus, there are superoscillatory incident wave functions which achieve the prescribed self- 
acceleration properties to arbitrary precision. For illustration, let us explicitly calculate such 
superoscillatory incident wave functions. 

5.1. Method of Matching Amplitudes. Let us begin by applying the method presented 
in jE], which is a special case of our method of general linear constraints. In this special case, 
we require the momentum-limited incident wave ip^(x) to exactly match the amplitude of 
the ideal template function at several points x k of the slit interval [ — §, §] • The constraints 
in the variational problem are then given by the linearly independent constraint functions 
Xkip) := and constraint parameters a k : = $(£&). Thus, 



1 



- far 



2nh 



m ax 



e f (** a* ) dp 



which leads to the solution: 



— — sine {— — [x k - x r \ ) (35) 



1 - 

^) {x)= J2\kXk(x), (36) 
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where A = T 1 a and where 

Xk(x) 



Pn 




e ft 



smc 



{x-x k ) 



dp 



Pu 



X 



(37) 



x 



is a linear combination of sine functions centred 



We observe that the wave function ip^ 
at the Xk and we note that ip(x) is square integrable, since the sine functions are. In general, 
T ill-conditioned, i.e. care must be taken to invert it with enough numerical precision so as 
to satisfy the constraints with sufficient accuracy. 

We used routines in Maple which calculate A = T _1 a by solving TA = a using Gaussian 
elimination. Concretely, we required ip( w '(x) to match the ideal template function <&(x) with 
p = 2 at N = 9 equidistantly-spaced points Xk from slit boundary to slit boundary. For 
example, FigH shows the imaginary part of the superoscillatory function ip( u > over the slit 

displaying the typical big amplitudes 



(«), 



x 



interval. FigEl shows a zoomed-out view of 
to the left and right of the superoscillating stretch. 



0.4 

/ \ 0,2 




1-3 j -2 \ -1 u l 
\s \ -0/2- 

\ / " 


i \^ 3 



Figure 4. Im (^^(x)) over the slit in the example of Sec l5.ll The wave- 
length is about 0.5A min . 

The momentum expectation value for the ideal template function that we chose is p = 
2pmax- Numerically, we found that the strictly momentum limited incident wave function 
ip(x) for N = 9 yields an emerging wave function ^(x) whose momentum expectation value 
is p ~ 1.92p max . Clearly, the momentum of particles which pass through the slit essentially 
doubles by self-acceleration, as intended. The momentum uncertainty of the emerging wave 
function is Ap lA2p max . 

Recall that for this slit size the momentum uncertainty could be significantly smaller, 
namely Ap = 1/2, as is precisely realized in the ideal template function. By increasing 
N, we can achieve that the incident wave function tp^(x) better matches the template, 
leading to a lowering of Ap towards that limiting value. For example, for = 15 we 
find p ~ 1.99947p max and Ap m 0.50025p max . For significantly larger A^ the exponential 
computational expense sets in. Our generalized method for linear constraints allows us to 
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Figure 5. ^(^(x)! over the slit and surrounding regions in the example of Sec l5.11 

use other linear constraints which we found to be numerically more efficient in the sense 
of allowing us to reach larger values of N. We will discuss the use of these alternative 
constraints in Sec l5.2l 




Figure 6. \^(x) — $(x)| 2 over the slit in the example of Sec 15. II 

FiglHl shows the accuracy with which the ip( u \x) obtained by matching N = 9 amplitudes 
of ip( u \x) to those of $(x) agrees with the ideal template function for arbitrary x in 
the slit interval. The behavior is similar for all values of N (that we tested). For general N, 
there are N — 1 peaks and the height of the highest peak decreases with N. 

5.2. Method of Matching Derivatives. In order to illustrate the generality of our new 
method of SecOl let us now construct superoscillatory wave functions by requiring that the 
wave function matches value and derivatives of the template at one point, instead of requiring, 
as we did in Sec l5.ll that the wave function matches only the value of the template function 
at several points. 

Concretely, let us require that the value and first N — 1 derivatives of the to-be-found 
wave function ip^ agree at x = with those of the 'ideal' template function $ of above. 
In the equation for general linear constraints, (|21j) . we obtain a constraint on the (k — l)st 
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derivative by choosing for the constraint function: 



Xk{p) 



ip 

T 



fc-i 



(3? 



Matching the derivatives to those of the template is to choose the constraint parameters to 
be ak := $^'~ 1 - ) (0), where <3>( fc_1 ) denotes the (k — l)st derivative. Since the Xk are linearly 
independent, 

1 /"P"" / ip^ -1 ( ip x r ' 



Ti 



kr 



2nh 



pmax 



dp 



(39) 



is invertible, yielding the solution 

^ {u \x) 



1 N 



fe=i 



where 



Pmax 



ip 

"h 



fe-i 



(40) 



(41) 



v2tt^ </-p max 

Note that ^(x) is a linear combination of derivatives of sine functions, each of which is 
bandlimited. In this case, T is simpler to invert and we can go, for example, to the case 
N = 23 before the exponential computational expense sets in. In this case, for large N the 
coefficients quickly grow large and hence the subtle cancellations in the Fourier transform 
require fast increasing numerical precision. 




Figure 7. \^(x) - $(x)| 2 over the slit in the example of Sec 15. 2 1 



We considered the example where the value of the function and its first 22 derivatives 
is required to match those of the ideal template function at x — 0. We found numerically 
that the momentum-limited superoscillating function, after passing through the slit, then 
exhibits a momentum expectation value of p ~ 2.0002p max and momentum uncertainty 
Ap fa 0.50049. Thus we reach the targeted momentum-doubling self-acceleration, with a 
momentum uncertainty which is only marginally above the uncertainty relation limit Ap = 
1/2 for this slit size. Fig0 displays the accuracy \^(x) — $(x)| 2 . 
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6. A Momentum Space Method 

In position space, superoscillatory wave functions if)(x) generally possess a characteristic 
shape: rapid but small oscillations in the superoscillating stretch and a few large long- 
wavelength amplitudes shortly before and after. Do these states also possess a characteristic 
shape in momentum space? 

Let us consider, for example, the superoscillations obtained by prescribing oscillating am- 
plitude values afc at close-by points x^. We found that, in momentum space, such a state is 
a linear combination of plane waves exp(— ix^p): 



V>(p) :-- 



1 \ 

y=^ r=1 X r e H 

0, if \p\ > p D 



if \p\ < Pv 



(42) 



It appears, see e.g. FigjHl that these ip(p) generally possess small amplitudes in most of the 
momentum interval [— j9 m ax,Pmax], except for the near the boundaries ±p max . We calculated 
the Fourier transforms of a number of superoscillatory wave functions and observed this as 
a general feature. 

Thus, in momentum space, these superoscillations appear to be a linear combination of 
plane waves whose interference is close to being as strong as it can be, with the effect that 
the resulting function is of minimized norm. 




Figure 8. Re (V>(p)) of the incident superoscillatory wave in the example in Sec l5.11 

If this assumption is correct, we should be able to derive superoscillatory wave functions 
by calculating that linear combination of plane waves in momentum space whose norm is 
minimal. To this end, let {x r }^ =1 be points in [— ■§,§]• Our aim is to find a coefficient 
vector {q r }r=\ of fixed length, say = 1 such that 



if \p\ > p n 



0. 



if \p\< Pi, 



(43) 



is of minimum norm. The constrained optimization problem with Lagrange multiplier v 



C 



+ Klk1 2 -i) 



(44) 
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leads to: 



Tq = vq 



(45) 



Thus, the coefficient vector q which solves this optimization problem must be eigenvector to 
T. From EqlH()| we obtain the general expression for the norm: ||?/>|| = q^Tq. Thus, q must 
be that eigenvector of T with the smallest eigenvalue. 

Indeed, the position wave function determined by these coefficients q is superoscillatory: 
already in Sec J4.2l we encountered the wave functions whose coefficient vectors A are the 
eigenvectors of T~ l of largest eigenvalue. There, we found that these are the superoscillatory 
wave functions which for a given set of constraint points {xk} are most norm-expensive and 
which, therefore, possesses the most pronounced superoscillations. 



We know from Sec J2.3.3l that it is always possible to find incident wave functions of fixed 
momentum bound that in the slit interval are matching any given template function arbi- 
trarily closely in the L 2 norm topology. Thus, for all practical purposes, the self-acceleration 
phenomenon can be tailored at will. Our method of general linear constraints can be used to 
explicitly construct a sequence of momentum-limited superoscillatory wave functions iPn(%) 
which more and more closely match any given template function $. The iPn{%) approach 
in the slit interval in the sense that they obey more and more linear constraints that 
tie iPn{ x ) t° 

In Sec l5.ll and Sec J5.2l we showed that a close approach to a fixed template function can 
be done numerically efficiently. Clearly, intuition and the easily achieved numerical accuracy 
lead us to conjecture that our methods for producing superoscillations, as used in Secs l5.1l5~2l 
do indeed always lead to convergence in the L 2 -topology towards the template function. So 
far, however, we have no proof that our particular method for producing superoscillatory 
wave functions from linear constraints does indeed realize the L 2 convergence to generic 
template functions. 

7.1. Quadratic constraints. Let us ask, therefore, if there is a choice of linear constraints 
that directly targets the area under the functions and that thereby directly guarantees con- 
vergence in the L 2 sense. 

One may try, for example, constraints which require that the functions if)( u ' and $ enclose 
equal areas on certain subintervals of the slit. This can be put into the form of a linear 
constraint: Let {x^}^ 1 be equidistantly-spaced points in [— §■>§]■ We require the linear 
constraints of (|21|) with the constraint functions 



7. Open Problems 




(46) 



and the constraint parameters: 




(47) 
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While this can easily be carried out, these constraints are not directly guaranteeing L 2 
convergence towards the template function by refining the partition of the slit interval into 
increasingly smaller subintervals: in principle, even functions that enclose equal areas on a 
very small interval may have very different amplitudes. Let us, therefore, consider to impose 
constraints which require the area of the function \^(x) — $(x)| 2 on small subintervals to 
be small. It is clear that to this end it would be necessary to implement also constraints 
that are quadratic in the field We will here not pursue this strategy to the end. As 

a preliminary step, however, let us generalize our method for constructing superoscillatory 
functions to include quadratic constraints. 

To this end, we formulate the variational problem of finding the function ip^ of smallest 
norm and with momentum cutoff p max which satisfies N linear and M quadratic constraints 
that tie it to a template function $: 



V2nh 
1 



V2nE 



Pmax 

^ {U) (P) XI(P) dp for k = l,...,N (48) 

Pmax 
Pmax 

ip {u) *{p) *p (u) (p) E* k (p) dp for k = 1, M (49) 



max 



The to-be-optimized functional with Lagrange multipliers {\k\h=\ and {^k}k=v rea ds 




L = I V (W) *(P) ^ {U) (P) dp~J2 -7= / ^ U) (P) Xl(p) dp 

i) {u) \p) i) {u \p) Et(p) dp + c.c. (50) 



in a x 



and the Euler-Lagrange equation reduces to 

- yfa z-mp) (bi) 

Although this may be difficult in practice, in principle, the substitution of (|51jl into (|4*H|) 
and (jUl)) yields sufficient equations to solve for the {\k)k=i and {fik}if=i in terms of the 
{ a k}k=i an d {^fcjfcLi and this yields the solution ip( u \ 

7.2. A Conjecture. Consider the case of a differentiable template function $ whose de- 
rivative is bounded: ^'(^Ol < ^Vi G [— fjfl, for some finite K. Assume that ip^ is 
a sequence of incident wave functions, calculated through the method of Sec 15.11 with the 
amplitudes of ip^\x) and $(x) matched at N equidistantly-spaced points x^. We conjecture 
then that the supremum ^^(x)! for all x and all N is finite as well: 



W N {x)\<M Vie 



L L 
2 1 2 



(52) 
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for some finite M. This is plausible because, else, would have to develop arbitrarily 

sharp spikes away from the template function in between some two points where its ampli- 
tudes are matched to those of the template function. From Sec l4.3| however, we expect large 
oscillations in the superoscillating stretch to be norm-expensive and therefore be prevented 
from occurring, given that the are optimized to possess minimum norm for a given set 
of constraints. 

7.3. Convergence. 

Proposition 2. Assume that the conjecture of Sec \7~^ holds true. Then, {iPn{x)}n converges 
uniformly and in the L 2 topology over the interval [— -|, \ ] to for N — > oo. 

Proof. Partition the slit into (N — 1) equal- length intervals with the N endpoints x^ : = 
-- + (&- 1)^. Define {a N (x)}% =2 by 

«iv(x) := max{ x^ \ k 6 {1, . . . ,N }, x^ < x }. (53) 

That is, a^(x) is the closest point in the partition from the left to x. Then, 

\ip N (x) - $(x)| < \ip N (x) - <$>(a N (x))\ + \$(a N (x)) - $(x)| 
= \i> N (x) - tp N {a N {x))\ + \<5>{a N {x)) - 
< M \x — ajv(x)| + K \o.n{x) — x\ 

<- ^ (M 

where we applied the triangle inequality and the mean value theorem. We therefore have 
uniform and L 2 -convergence. □ 



8. Summary 

We started with the method for calculating superoscillatory wave functions introduced in 
[U] and applied it to concrete examples. We then generalized this method so that it now 
allows us to construct superoscillatory low-momentum wave functions with a wide range 
of predetermined properties. Namely, we can impose any arbitrary finite number of linear 
constraints. We calculated concrete examples. 

Further, we addressed the question whether superoscillatory functions can be made to 
match any arbitrary continuous function on a finite interval. This would correspond to 
imposing an infinite number of constraints. Generally, the answer is no. However, we were 
able to prove that there always exists a sequence of superoscillatory wave functions which 
converges in the physically relevant L? topology towards any continuous template function 
over an arbitrarily large chosen interval. 

This is of interest for example in the case of the single slit: we proved that the wave 
function of an incident low-momentum particle can be chosen to arbitrary precision such 
that, if the particle passes through the slit, it will emerge with a predetermined arbitrarily 
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large momentum expectation and with a momentum uncertainty that is as small as permitted 
by the width of the slit. 

Our method for constructing superoscillating wave functions allows us to construct super- 
oscillatory wave functions which match any finite number N of properties of a given template 
function (such as the template function's amplitudes or derivatives at specified points). This 
leads to the question if by letting the number of constraints, N, go to infinity we can ob- 
tain one of those sequences of superoscillatory wave functions which converge towards the 
template function in the L 2 topology. 

We proved that such sequences exist but we have not proved that our particular method 
produces such sequences. The numerical evidence certainly suggests that this is the case. In 
fact, we found rather fast numerical convergence. 

Nevertheless, it would be highly desirable to be able to prove that a given method for 
producing superoscillations can be used to calculate a sequence of superoscillatory functions 
that converges in the L 2 topology towards any given template function on an interval. An 
investigation based on Weierstrass' approximation theorem is in progress, [T5] . 

Lastly, we found a method for identifying a class of superoscillatory functions by looking 
at their behavior in momentum space: superoscillatory functions can be viewed as functions 
which in momentum space are a linear combination of plane waves with coefficients such 
that their interference is maximal, i.e. such that their norm is minimal. 
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